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(N ; ABSTRACT 

C3 t An analysis of cosmological mixed dark matter models in spatially flat Friedmann Universe with zero A- 
term is presented. We argue that the introduction of cosmic gravity waves helps to satisfy observational 
constraints. The number of parameters is equal to 5, they are (1) as, the dispersion of the mass fluctuations 
in the sphere with radius 8 h~ 1 Mpc, (2) n, the slope of the density perturbation spectrum, (3) Q u , the density 
£N ■ of hot dark matter, (4) fib, the density of baryons, and (5) h, the Hubble constant Hq = 100 h kms^Mpc^ 1 . 
The cold dark matter density parameter is equal to Sl c d m = 1 — &v — The analysis of models is based on 
the confrontation with the mass function of clusters of galaxies and the CMB anisotropy. The implication 
| of Press-Schechter formalism allowed to constrain a% = 0.52 ± 0.01. This normalisation of the spectrum of 

■ density perturbations has been used to calculate numerically the value of the large scale CMB anisotropy and 
. the relative contribution of cosmological gravitational waves, T/S. We found that increasing £l u weakens 
| the requirements to the value of T/S, however even for £2 U < 0.4 the models with h + n > 1.5 suggest 

■ considerable abundance of gravitational waves, T/S^0.3. In models with Q u < 0.4 and scale-invariant 
^ spectrum of density perturbations (n = 1), T/S^10(/i — 0.47). Minimisation of the value T/S is possible 
O | only in the range of the red spectra (n < 1) and small h (< 0.6). However the parameter Vt u is strongly 

' constrained by AT/T data on the first acoustic peak of Sakharov oscillations. Assuming that T/S G [0,3] 
03 ■ and taking into account observational data on the primordial nucleosynthesis and the amplitude of the first 
acoustic peak we constrain the model parameters. We show that the considered models admit both moderate 
"k^J | red and blue spectra of density perturbations n G [0.9,1.2] with rather high abundance hot dark matter, 
j_j ■ Q u G [0.2,0.4]. Any condition, n < 0.9 or Q, u < 0.2, decreases the relative amplitude of the first acoustic 
. 5^ , peak for more than 30% in comparison with its height measured by BOOMERanG and MAXIMA-1. 



INTRODUCTION 

A great success of observational cosmology of the last decade has allowed to go from a purely theoretical 
investigation of possible cosmological models to testing them and constraining their parameters. Currently, 
the total number of model parameters exceeds ten, they can be separated in few groups: 

• The parameters related to the geometry of the Universe - QtotaU Hubble parameter Hq, and three- 
dimensional curvature k. 

• The parameters related to the nature and abundance of matter (including dark matter and baryons) 
in the Universe and the value of cosmological constant (A-term), called often as "the dark energy". 
They are the matter density parameters in the Universe ^matter (the sum of £l c dm, and f^), the 
number of massive species of neutrinos N v , the energy density of A-term Q\, and other parameters of 
the dark energy (e.g. quintessence, its equation of state p = wp, etc). 

• The group of parameters related to the inflationary stage of the Universe, which describes the power 
spectra of metric perturbations generated. In the simplest version of inflation based on the theory of 
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scalar field these parameters are the amplitudes and the slopes of density perturbation and gravitational 
wave spectra. 

A complete analysis of the cosmological model in constraining all mentioned above parameters (assumed 
initially as free ones) is currently a difficult task since it requires as excellent technical facilities as a special 
study of different sets of observational data for their completeness and independence. In view of the im- 
perfection of currently available observational data one usually assumes a simplified model with restricted 
number of the parameters (and fixed others) and test it by a selected set of observational data. 

The goal of the present analysis is to consider a matter-dominated model of the Universe and constrain its 
parameters by three fundamental tests which have different characteristic scales. These observational tests 
are the large scale CMB anisotropy (~ 1000 h~ 1 Mpc), the amplitude of the first acoustic peak of Sakharov 
oscillations (~ 100 h~ 1 Mpc) and the mass function of clusters of galaxies (~ 10 h~ 1 Mpc). The aim of the 
analysis is to clarify the modern status of such a model: whether it can explain the observations or not, and 
thus, whether the ACDM is the only satisfactory model. 

To answer this question we considered a family of matter-dominated cosmological models with the fol- 
lowing free parameters: 

• cr 8 G [0.47,0.61], (15 models with step 0.01); 

• n £ [0.8, 1.4], (7 models with step 0.1); 

• fij, € [0,0.4], (5 models with step 0.1); 

• Q b G [0.01,0.11], (6 models with step 0.02); 

• he [0.45,0.7], (6 models with step 0.05). 

Altogether we had 18900 variants of the model. The broad intervals for parameter variations are specially 
selected to find all possible solutions for the given tests. 

The first parameter, as, is a more accurate normalisation of density perturbation spectrum than COBE 
one, since a probable contribution of primordial gravitational waves into the COBE measured signal. The 
latter is the only available today data which we can use to find a signature of cosmic gravitational waves. 

The second parameter, n, generalises the scale invariant spectrum of density perturbations (Harrison- 
Zeldovich spectrum) in the simplest way. A power-law spectrum with n < 1 (the "red" one) is predicted in 
a broad range of inflationary models (Linde, 1983, Lucchin and Matarrese, 1985). In other models it can 
be both "blue" (n > 1) and non-power-law (the noncomplete list of references includes: Carr and Gilbert 
(1994), Gottlober and Miicket (1993), Lesgourgues et al. (1997), Lukash and Mikheeva (2000), Lukash et 
al. (2000), Melchiorri et al. (1999), Mollerach et al. (1994), Semig and Midler (1996)). Here we consider 
the power-law spectra with a free slope index n within the range. 

The introduction of hot dark matter is related with a recent discovery of the atmospheric neutrino oscil- 
lations, which demonstrates that at least one species of neutrino is massive. A careful theoretical analysis 
of all neutrino experiments argues for three massive species, but this result still requires a confirmation. In 
the numerical analysis we assumed that N u = 1. Actually, the value of N v becomes important only at small 
scales (<1 h~ 1 Mpc, e.g. in the analysis of Ly a -absorbers, the test which we do not consider here). 

The fourth parameter, f^, influences the transfer function slightly, so, neither mass function of clusters 
of galaxies nor large scale CMB anisotropy can reveal its value with a high confidence level. Hovewer, Ctf, is 
extremely important for the amplitude of the first acoustic peak of Sakharov oscillations. 

The last parameter is Hubble constant. Till now there is some discrepancy in its measurements (Saha et 
al., 1999, Jha et al., 1999). It is clear that h cannot be smaller than 0.5, but on the other hand, the large 
value (>0.65) will contradict the age of the oldest globular clusters. 

In our analysis there are two derived parameters. The first of them is the abundance of cold dark matter, 
^■cdm = 1 — fij/ — The second one is a relative contribution of primordial gravitational waves into large 
scale CMB anisotropy, T/S. The latter parameter can be theoretically calculated in any inflationary model. 
T/S is not discriminated by inflation and vary in a broad interval (we assume T/Sg [0,3] to take into 
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Fig. 1. 7V(> M) as a function of a 8 (n = 1, Sl v = 0, 
£lb = 0.05, h = 0.5). The upper curve corresponds 
to the model with a% = 0.61, the middle one - as = 
0.52, the lower curve - as = 0.47. 



account all possible values). For the majority of inflationary models there is a simple relation between the 
spectrum slope of primordial gravitational waves, tit, and T/S: 



T/S ~ -6n T , 



(1) 



where the numerical coefficient slightly depends on scale ( T/S is calculated here for 10° angular scale). 
Eq.(l) allows to replace one parameter (T/S) by another («t)- It is interesting to note that in power- 
law inflation there exists an additional relation between n and tit, so in this model T/S can be expressed 
in terms of n. However, generally such relation is absent (see, e.g., Lukash et al. (2000), Lukash and 
Mikheeva(2000)). Thus, instead of considering additional constraints related with given inflationary model, 
we prefer to consider n and T/S as independent variables within their ranges. 

MASS FUNCTION OF GALAXY CLUSTERS 

The number of massive halos with the mass larger than M can be calculated with help of Press-Schechter 
(P&S) formalism (Press and Schechter, 1974): 



N(> M) 



dn 
dM 



M 



where 

dn J~2 pod, 
dM = " 



7 dM', 



da R 



tt M al 



R 



dM 



cxp 



(2) 



(3) 



where M = 4poR 3 /3, po is the matter density in the Universe, 8 C is a critical density contrast necessary to 
form a halo (5 C = 1.686, see Gunn and Gott (1972)), an is a dispersion of linear density perturbations in 
the sphere of radius R which is an integral of the power spectrum of density perturbations 



AiW'(kR)- 
o K 



(4) 



where A? 



i 



\ k — 2^zP(k)T (k) is a dimensionless power spectrum, P(k) = Ak n+3 - a primordial spectrum of 
density perturbations, A - a normalisation constant, T(k) - a transfer function which depends on 17^, 
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log 10 M, [h 'MJ 
Fig. 2. AT(> M) as a function of Vl v (a 8 = 0.52, n = 
1, $7b = 0.05, h = 0.5). The solid line corresponds to 
the model with Q u = 0, the dashed one - Vt v = 0.2, 
the lower line - Q u = 0.4. 



and h (we used the analytic approximation of the transfer function given by Novosyadlyj et al. (1999)), 
W(kR) - Fourrier transform of top-hat filter, W(x) = 3 (sinx — xcosx) /x 3 . 

The Figure 1 illustrates how the function N(> M) depends on <7g with other parameters fixed. As it was 
mentioned in the Introduction, erg affects efficiently the height of P&S curves which one can use to determine 
<7g. The Figure 2 demonstrates the behaviour of N(> M) as a function of Q, v . Taking a given o~% (0.52 in 
the Figure 2) fixes the amplitude of P&S curve at the mass value in the sphere with radius 8 h~ 1 Mpc, while 
different values of fi„ change the slope of the curve. In both Figures the circles with errorbars mark the 
observational data. For the rest parameters, N(> M) varies with a much smaller amplitude and we do not 
adduce the corresponding figures. The observational data in the Figures 1 and 2 are taken from Bahcall 
and Cen (1993). 

To analyse the cosmological models we have used the standard \ 2 method. At the first step we have 
calculated x 2 f° r eac h model (this function is x 2 distributed with 5 degrees of freedom; Nd.f. = N b s ,p i nts — 
^parameters = 10 — 5). It turns out that the model with the minimal x 2 (Xmin = 1-43) has the following 
parameters: as = 0.52, n = 1.3, Q, v = 0.3, f2& = 0.01, h = 0.7. To evaluate the confidence level of this result 
we took into account the mathematical theorems available for x 2 an d constructed other distributions for 
each parameter A%f (a) = x\~ X m in> where a is one of the five cosmological parameters, x\ i s a minimal 
X 2 on the hypersurface a = const in 5-dimensional model space. Ax 2 (a) is also x 2 distributed but has the 
only degree of freedom. The confidence levels 1, 2 and 3a correspond to Ax 2 (a) = 1, 4, and 9 respectively. 
Ax 2 (a) for a = as is shown in the Figure 3 (the curve marked by " x "). One can see that as = 0.52 ± 0.01. 
The uncertainties of the Press-Schechter formalism and data systematics enhances the total errorbar to 0.04 
(Eke et al., 1996, Borgani et al., 1999). 

To evaluate the accuracy of the approximation of the transfer function and its influence on the output 
we have repeated this statistical analysis (see Figure 3, the line marked by " + ") with another analytic 
expression for the transfer function (Eisenstein and Hu, 1999) and found that both approximations lead to 
the same result. 

As to the other cosmological parameters (n, £l u , h) their Axl(a) are under la confidence level, 
therefore they cannot be constrained by the cluster data only. So, as a result of this part of our analysis 
we have obtained the high accuracy normalisation of the density perturbation spectra. It is evident that the 
obtained value of as depends on observational data used and can be slightly different if one uses other data 
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The curve marked by "+" corresponds to the calcula- 
tions with the transfer function of Eisenstein and Hu 
(1999). The curve marked by " x" corresponds to the 
transfer function from Novosyadlyj at al. (1999). The 
horizontal lines mark the different confidence levels. 



(e.g. Borgani et., 1999). 

CMB ANISOTROPY 

Large Scale CMB Anisotropy 

The large scale CMB anisotropy is related with metric perturbations (Sachs and Wolfe, 1967): 



AT.. 1 f E dh ik 



e i e h dr], 



(5) 



where E and R are emission and reception moments of time, is a metric perturbation tensor, d/drj is 
the conformal time derivative, the integral goes along the null geodesies. 

The relative contribution of cosmic gravitational waves into the large scale CMB anisotropy (T/S) can 
estimated as follows: 



S+T=S 1+ 



10° 



~ 1.1 x 10 



-10 



(6) 



where S is the contribution of the density perturbations: 



S = J2 SeWi, 



1=2 

where 



a _2l+l, ffra +3 r(3-n)r(l+(n-l)/2) 
1 64^ T 2 (2-n/2)T(£ + (5-n)/2)' 



and 

We = exp 



2^ + 1 
27 



(7) 



(8) 



(9) 
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4. T/S as a function of n and fl u (og = 0.52, 
0.05, h = 0.5). 



We evaluate the normalisation constant A from Eq.(4) for a 8 , Wi is COBE DMR window function. 

Figure 4 illustrates T/S as a function of n and Q. v for as = 0.52, Slf, = 0.05 and h = 0.5. The thick line 
represents the T/S-n relation for the power-law inflation. It can be seen that small values of T/S favour red 
spectra of density perturbations. 

The Figure 5 demonstrates T/S as a function of n and h for ag = 0.52, = 0.05, and large Q u = 0.4. 
The value T/S increases linearly with h and decreases with Q, y growing, therefore the curves T/S with the 
maximum parameter Q u = 0.4 can be used to put the lower limit on T/S (see Conclusions). 

Taking a moderate T/S< 0.5 and nearly flat power spectrum (0.92 < n < 1.02), we put an upper limit 
on the Hubble constant, h < 0.6, and lower limit on the hot dark matter abundance, fi„ > 0.1. However, 
the hardest constraint for the parameter Q. v can be got when we confront the amplitude of the first acoustic 
peak in AT/T (t ~ 200) with the observational data. 



For <7g = 0.52, SI;, = 0.05 and h = 0.5 we propose the following approximation for T/S (the accuracy is 
better than 11 % for 0.1 <T/S< 3): 



The Medium Scale CMB Anisotropy 

To compare the height of the first acoustic peak generated in our models with its value measured by 
BOOMERanG (Lange et al., 2000) and MAXIMA-1 (Balbi et al., 2000) we introduce the dimensionless pa- 
rameter 5ft which is the amplitude of the peak in comparison with the large scale CMB anisotropy amplitude: 
3? = 3fy=20o/l-l x 10~ 10 , where = 1(1 + l)S e /(£ + 0.5); 5ft ~ 5 for BOOMERanG. 

Evidently, 3? decreases with T/S growing (and other parameters fixed). E.g. for normalised by as standard 
Cold Dark Matter models (Q y = 0, n = 1, flf, = 0.08, h = 0.5) the relative amplitude of the acoustic peak 
decays by a factor T/S+l ~ 4 and practically disappears. 

The more efficient ways to enhance the first acoustic peak is a transition to the red power spectra and/or 
high abundance of the hot matter. The role of the blue spectra becomes important when the parameter S7„ 
rises up (since the red spectra will contradict to the condition T/S> 0). The results of the derivation are 
presented in the Figure 6. The condition for a "considerable" acoustic peak (5ft > 3.5) with the instantaneous 
Big Bang nucleosynthesis constraint for the baryonic density, leaves us with the standard family of the power 
spectra (n € [0.9, 1.2]) but requires high fraction of the hot matter (Q, v € [0.2, 0.4]) in the class of the models 
considered. 

Varying we found that increasing Q b by factor of two results in a significant rise of the acoustic peak 
independently on n and Q u , so for models with = 0.1, n £ [0.9, 1.2] and Q u € [0.2,0.4]: 4 < 5ft < 5.5. As 




(10) 
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to the parameter h, its increasing decreases the acoustic peak amplitude. Summarizing, we can conclude 
that MDM models with gravitational waves can marginally reproduce the BOOMERanG result indeed. 

CONCLUSIONS 

• The data on the galaxy cluster abundance determine the value as with a high accuracy, the other 
parameters (n, Q, u , h) remain free within their ranges. 

• None of the MDM models with n = 1 and T/S = satisfies both normalisations, on the galaxy cluster 
abundance and large scale AT/T anisotropy, which leads either to rejection from the flat spectrum or 
to the introduction of a non-zero T/S (or both). 

• Small values of T/S are realised for the red spectra (n < 1) and moderate h(< 0.6), the violation of 
these conditions leads to a high T/S(>1). 

• Increasing Q u weakens the requirement to the value of T/S, however even for f2„ < 0.4 the models 
with h + n > 1.5 suggest considerable abundance of gravitational waves: T/S^0.3. 

• In models with fl u < 0.4 and scale-invariant spectrum of density perturbations (n = 1): T/S>10(/i — 
0.47). 

• In models with = 0.05 and h = 0.5 we have the approximation formula for the primordial gravita- 
tional wave abundance (see Eq.(10)). 

• In double-normalised models with T/S>0 the height of the acoustic peak is less than its 'standard' 
value (3? = 5). The decrease of the parameter 3? does not exceed 30% in models with large fi„ 
(€ [0.2,0.4]) and any spectrum slope, n £ [0.9, 1.2]. Any condition, n < 0.9 or Q u < 0.2, decreases the 
relative amplitude of the first acoustic peak for more than 30% (i.e. 5? < 3.5 in models with = 0.05, 
h = 0.5). The acoustic peak practically disappears in CDM models. 

• When growing the baryonic abundance the factor 3? increases. The amplitude of the acoustic peak 
coincides with its 'standard' value (3? = 5) up to accuracy 10% in models with 0,b = 0.1 and both, 
blue spectra n G [1, 1.2] and Q u > 0.3 as well as moderate red spectra n <G [0.9, 1] and Q. v > 0.2. 

• Thus, rising the parameter £l u up to the values in the interval [0.2, 0.4] solves the problem of the first 
acoustic peak in AT/T, leaving the baryonic density within the primordial nucleosynthesis constraints. 
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Fig. 6. 5ft as a function of n for different values of 
n v {n b = 0.05, h = 0.5, o- 8 = 0.52). Models in the 
non-shaded region have T/Se [0,3]. 
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